Theory of Transient Phenomena in the Transport of Holes 
in an Excess Semiconductor 

By CONYERS HERRING 

An analysis is given of the transient behavior of the density of holes m, in an 
excess semiconductor as a function of time / and of position x with respect to the 
electrode from which they are being injected. When the geometry is one-dimen- 
sional, an exact solution for the function >th(x, t) can be constructed, provided 
certain simplifying assumptions are fulfilled, of which the most important are that 
there be no appreciable trapping of holes or electrons and that diffusion be negligi- 
ble. An attempt is made to estimate the range of conditions over which the 
neglect of diffusion will be justified. A few applications of the theory to possible 
experiments are discussed. 

A variety of experiments have been performed, and others are planned, 
which involve measurement of transient or steady-state phenomena due to 
the drift of positive holes along a specimen of //-type semiconductor after 
they have been introduced at an injection electrode or emitter. 1 These phe- 
nomena are presumably a result of the interplay of drift, space-charge, re- 
combination, and diffusion effects. This paper seeks to relate these effects to 
the phenomena, and its principal contribution is an explicit calculation of the 
transient phenomena outside the range of small-signal theory, for cases 
where the geometry is one-dimensional and where certain simplifying as- 
sumptions, notably the neglect of diffusion, are justified. Removal of some 
of these simplifying assumptions and a more careful development of the 
theory will be necessary in certain applications. 

Section 1 discusses the physical assumptions and boundary conditions 
involved in setting the problem up. Section 2 contains calculations of the 
distribution of holes along the length of the semiconductor at various times, 
for the mathematically simplest case where recombination and diffusion are 
ignored and all currents are held constant after the start of the injection. 
This simple case illustrates the method of attack, to be used in the more 
general calculations of Section 4, and it is hoped that this sketching of basic 
ideas will enable the hasty reader to pass on to Section 6 without going 

1 Experiments of this sort have been undertaken with the objective of testing and 
extending the theoretical interpretation of transistor action proposed by J. Bardeen and 
W. H. Brattain, Phys. Rev., 75, 1208 (1949), especially as regards the role of volume 
transport of holes, a role first suggested by J. N. Shive, Phys. Rev., 75, 689 (1949). Ex- 
amples of the type of experiment discussed in the present paper have been described by: 
J. R. Haynes and W. Shockley, Phys. Rev., 75, 691 (1949) (transient effects); W. Shockley, 
G. L. Pearson, M. Sparks and YV. H. Brattain, in a paper presented at the Cambridge 
Meeting of the American Physical Society, June 16-18, 1949 (steady-state transport); 
W. Shockley, G. L. Pearson, and J. R. Haynes. Bell Sys. Tech. Jour., this issue (steady- 
state and transient effects). 
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through the mathematical details of Sections 3, 4, and 5. Section 3 contains 
the complete differential equations of the problem, including diffusion and 
recombination, and Section 4 gives the solution when only the diffusion 
terms are neglected. Section 5 contains some order-of-magnitude estimates 
regarding diffusion effects. Section 6 summarizes the capabilities of the 
theory so far developed, presents some obvious generalizations, and dis- 
cusses an interesting shock wave phenomenon which occurs whenever the 
injected hole current is quickly decreased. 

1. Basic Assumptions and Boundary Conditions 

Consider the w-type semiconducting specimen shown in Fig. 1, having 
electrodes at its two ends, * = -a and :v = b, respectively, and an injection 
electrode system at * = somewhere in between. Let a current of density j a 
per unit area enter at the left-hand end, and let a current of density j e be 
injected at x = 0. To make the problem strictly one-dimensional, it will be 
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Fig. 1— Idealized experiment on hole transport in one dimension. 

supposed that this injection takes place uniformly over the plane cross- 
section of the specimen at x = 0, instead of taking place at isolated points 
of the surface, as is usually the case in experiments. This idealization will 
presumably be justified if the thickness of the specimen is small compared 
with lengths in the ^-direction which are significant in the experiment and 
if the injected positive holes are able to spread themselves uniformly over 
the cross-section before appreciable recombination has taken place. 

Unless otherwise stated, it will be supposed that j, consists entirely of 
positive holes, i.e., that the number of electrons withdrawn from the speci- 
men by the electrode at % = is negligible compared with the number of 
holes injected. The currents j a &ndj e need not be constant in time, although 
most of the analysis to be given below will assume them constant after the 
time of initiation of j e . 

One can set up differential equations for the variation with x and time 
of the electron density, n, , and the hole density, n h . These equations will in 
the general case involve migration due to electrostatic fields, diffusion, re- 
combination, trapping, and thermal release of electrons and holes from 
traps. It will be assumed, however, that trapping and thermal release from 
traps can be neglected, or, more precisely stated, that creation of mobile 
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holes and electrons occurs only at the electrodes, and that the disappearance 
of mobile holes and electrons is caused only by mechanisms which cause 
holes and electrons to disappear in equal numbers at essentially the same 
time and place. If this assumption is valid, the charge density due to im- 
purity centers will never differ from its equilibrium value by an amount 
comparable with the density due to free electrons. This assumption can be 
expected to be reasonably good for an w-type impurity semiconductor in 
which the number of donor levels is very much greater than the number of 
acceptor levels and for which, at the operating temperature, practically all 
the donor levels have been thermally ionized, while thermal excitation of 
electrons from the normally full band has not yet become appreciable. 

As has just been mentioned, the differential equations for the behavior 
of the electron and hole densities involve migration under the influence of 
the local electric field E(x, I). This field is in turn influenced by the space 
charge due to any inequality between the hole density n h and the electron 
excess (n e — » ), where n is the normal electron density. If the difference 
(»/, — n, + n ) were comparable with tin or n e , the problem would be very 
complicated. Fortunately, however, this difference cannot have an appreci- 
able value over an appreciable range of x, on the scale of typical experiments. 
For example, if {n h — n e -f- «o) were 10 -2 of n for a range Ax of 1/x, and if 
«o is 10 15 cm -3 , then the difference in field strength on the two sides of Ax 
would be about 2000 v/cm, a field which would outweigh all other fields in 
the problem and rapidly neutralize the space charge. Moreover, the time 
required for the evening out of any such abnormally high space charge would 
be very short, of the order of magnitude of the resistivity of the specimen 
expressed in absolute electrostatic units (1 sec. = 9 X 10" fi cm). Thus it 
will be quite legitimate to assume (tih — n e + n ) = in all equations of 
the problem except Poisson's equation which determines the field E, and 
so n e can be eliminated from the conduction-diffusion equations for holes 
and electrons. These two equations can then be used, as is shown below, to 
determine the two unknown functions tin and E, Poisson's equation being 
discarded as unnecessary. 

The boundary conditions for these differential equations consist of two 
parts, the conditions at I = and those at and to the left of x = 0. In most 
of the applications to be considered, the injection current j',. will be assumed 
to commence at / = 0. Thus, initially, the specimen will be free of holes and, 
at / = 0+ will have a field E a = j a /<r in the region — a < x < 0, and a 
field E = jb/<To in the region < x < b, where a is the normal conductivity 
of the specimen and j b = j a + j e is the total current density to the right 
of x = 0. The boundary condition at x — is determined by the magnitudes 
of the electronic and hole contributions to the injection current j e . If no 
electrons are withdrawn by the electrode at x = 0, then the electron cur- 
rents just to the left and just to the right of x = must be equal, and the 
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hole current densities on the two sides must differ by j ; if a part of j e is 
due to withdrawal of electrons, then the electronic current will have a cor- 
responding discontinuity. Uj a is positive, i.e., flows from left to right in the 
specimen, the current can be assumed to be practically entirely electronic 
over most of the range from -a to 0; i.e., as x becomes negative the hole 
current must rapidly approach zero and the electron current must rapidly 
approach j a . In fact, if diffusion is ignored the electron and hole currents 
must have these limiting values for any negative x. 

The preceding discussion and the mathematics to follow have been 
couched in purely one-dimensional language, i.e., have been formulated as if 
the electron and hole densities were functions of * alone, independent of y 
and z, and as if the semiconductor extended to infinity in the y- and z-direc- 
tions. However, it is easy to see at each stage that practically the same 
equations can be written for transport of holes along a narrow filament whose 
thickness is small compared with the linear scale of the phenomena along its 
length, even when the density of holes' is not uniform over the cross-section 
of the filament. If the density of holes is uniform over the cross-section, all 
the equations will of course hold as written. However, recent work 2 has 
suggested that holes recombine with electrons so rapidly at the surface that 
the density of holes may be much smaller near the surface than in the center 
of the cross-section. In such case all the equations of this memorandum must 
be interpreted as applying to the mean value, ii h (x), of the density of holes, 
nh (x, y, z), averaged over the cross-section of the filament; also, the rate of 
recombination of holes and electrons must be set equal to some function of 
n h , as yet not reliably known, instead of to a constant times the product of 
electron and hole densities. This will of course alter most of the quantitative 
predictions of Section 4, but will not require any change in the method of 
calculation. 

2. Formulation and Solution of the Problem with Neglect of 
Diffusion and Recombination 

For this case the electron and hole currents can each be equated to the 
product of field strength E by particle density n by mobility m, and the 
continuity equations are 

^ =-»(/ 4tt ,„) (0 

at ox 

a 3=i (£„,„,). (2) 

dt dx 

• H Suhl and \V. Shockley, paper Oil presented at the Washington Meeting of the 
American Physical Society, April 29, 1949; see also Shockley, Pearson, Sparks and Brat- 
tain, reference 1. 
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Since the neutrality condition requires — = — -, subtracting(l)and(2)and 

ot at 

integrating gives the equation of conservation of total current: 

E(jt e n e + MftWft) — j(0/ e = const, indep. of x 

where of course j = jb = (j„ + /«) when < x < b and when conditions 
are such that all currents flow from left to right. Putting the neutrality con- 
dition n e = )ih + »o , into the equation gives the following relation between 
E and «/, : 

E[(n c + nh)n h + fieHo) = j/e (3) 

This can be used to eliminate either E or »* from (1). If £ is eliminated 
we have 

dl e[(n e + Ma)"/. + Me«o] 2 dx dx 

where V(ii-h) is an abbreviation for the coefficient shown. If, instead, Hi, is 
eliminated from (1) a similar equation results: 

«-||- F «)« (5) 

dl j at ox 

where 

V(E) = cE-y-kHcUa/j = Ep h (E/Eo) (6) 

where 

E = j/<rt (7) 

i.e., the field necessary to maintain the total current by electronic conduction 
in the normal state of the specimen. The velocity V(E) is of course numeri- 
cally the same as the V(uh) occurring in (4) when E and n h are related 
by (3). 

The solution can be based on either (4) or (5). We shall use (4), as »/, 
is the most interesting quantity for direct measurement, and as the differen- 
tial equation to be given below for the case where diffusion terms are in- 
cluded is simpler when nn is chosen as the dependent variable. 

Equation (4) (or (5)) describes a wave propagated with the variable 
velocity V. If j e « j a , so that E is never greatly different from E , 
(4) (or (5)) and (6) indicate that m (or E) is propagated with the constant 
velocity E ma > as is of course to be expected. More interesting is the case 
where /, and j a are comparable, so that V departs significantly from con- 
stancy. It is tempting to suppose that, for this case also, the curve of tin 
against x at any time / can be constructed by taking the graph of iik against 
x at / = and moving each point of the curve horizontally to the right a 
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distance V(n h )t. One can, in fact, easily verify that this construction gives a 
solution of (4), by writing (4) in the form 

(dni\ 

\dt/n h fdnk\ 

\dx)t 

whence it is obvious that the function n h {x, t) defined implicitly by 

x(n h , /) = x(n h , 0) + V(n h )t 

satisfies (4) for any form of the arbitrary function x(n h , 0), and that, con- 
versely, any solution of (4) must be of this form. 

Assuming, as in the preceding, that all currents flow from left to right, 
the boundary conditions at / = + are : 

n h = for x < and x > (8) 

or, equivalcntly, 

E = E a = j a /c for * < Oj 

E = E - (Ja + /•)/«* for x > Oj 

The boundary conditions at x = are, for t > 0, 

n h = or, equivalently, E = E a for X = 0" (10) 

and 

»a = »Ai or, equivalently, E = Ei , for x = 0+ (11) 

where £i and n h \ are given by the requirement of continuity of electronic 
current, i.e., 

Ealing = Ei(n + Uhljfic 

whence, using the relation (3) between £i and » u and expressing £ a as 
j a /noey. t 



n h i = 



n ° (12) 



or, alternatively, 



■Z?^f* _ i 



a.jfcfi-fa+j^ *^1 (13) 

According to (12), » M is small when j e is small; and, by (13), £i is only 
slightly below £ for this case. As/ e increases, n h i increases and E x decreases, 
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and (12) and (13) would make »« infinite and E% zero when j e /j a .= aVV« » 
This merely means that the assumptions made in this section, in particular 
the neglect of diffusion and recombination or the assumption that no. elec- 
trons are taken out by the injection electrode, must fail to be valid before 
j c gets as large as ma>/mo • It will, in fact, be shown in Section 5 how the 
presence of enormous concentration gradients makes it essential to consider 
the effects of diffusion near x = when j e becomes large. 
Putting the boundary conditions (8), (9), (10), and (11) into the wave- 
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Fig. 2 — Schematic variation of hole density n\ and electric field E with distance x from 
injection electrode and time / after the start of the injected current, in the approximation 
neglecting diffusion and recombination. 



propagation construction described above gives the solution shown schemati- 
cally in Fig. 2. An infinitesimal instant after / = 0, rtj, is zero everywhere 
except in an infinitesimal interval at x = 0, where it rises to a maximum 
value iihi given by (12). This is shown schematically in the upper left dia- 
gram of Fig. 2. The corresponding plot of E, shown in the upper right, dips 
down to Ei , which is less than either E a or E , in this infinitesimal interval. 
After a finite time has elapsed, the curves of ii/.i and E against x are simply 
those obtained by moving each point of the right-hand portions of these 
/ = + curves a distance VI horizontally to the right, as shown in the bottom 
two sketches. Here V depends on the ordinate in each diagram, taking on 
its maximum value E ui t when Vh\ = or E = E . Since V is proportional 
to E-, the curve in the lower right diagram is a parabola in the range be- 
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tween the front and the rear of the transient disturbance; this parabola, if 
continued, would have its vertex at the origin. After a sufficiently long time 
a steady state will be reached in which the field for positive x has the uni- 
form value Ei and the density of holes the uniform value iihi . 

It is possible to measure n h as a function of / for fixed x by using a closely 
spaced pair of probes to measure the potential gradient E, and converting E 
to n h by (3) ; alternatively, the current to a single negatively biased probe 
can be used as a measure of n h , if calibrated by the two-probe method. The 



nhi 




Fig. 3— Schematic variation of hole density if* with time / after the start of the injected 
current, at some given distance downstream from the injection electrode, in the approxi- 
mation neglecting diffusion and recombination. 

portion of this curve of n h against / for which < n h < n h i is given, in the 
present approximation, by 



where 



= / F [l + (1 + nu/tie)>ih/no]- 



//, = x/Eon, t 



(14) 



(15) 



is the time of arrival of the front of the disturbance. This curve is a parabola, 
as shown in Fig. 3; if continued, the parabola would have its vertex on the 
negative n h axis, as shown. The rear of the disturbance, at which n h becomes 
constant and equal to n h \ , arrives at a time t R given by inserting n hl from 
(12) into (14): 

i R = /,/u - (i + nMWUa + y.)i 2 ( 16 ) 
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Note that the velocity of advance of the rear of the disturbance is less than 
that with which the holes drift in the steady-state field E x . In other words, 
wave velocity and particle velocity must be distinguished in phenomena of 
this sort, although they happen to coincide at the front of the disturbance. 

The discussion just given has been based on the assumption thatja and 
/Vare independent of time, and that they both flow from left to right in Fig. 1. 
Time changes in the currents are easily taken into account in the iii, con- 
struction of Fig. 2: according to (4), it is merely necessary to move the 
various points of the curve of tin against % to the right with the variable 
velocity V(tth , instead of the constant velocity V()ih); in addition, m,i 
will in general not be a constant, so that the part of the curve for small x 
will no longer be a horizontal line. As for the restriction that the currents 
all flow from left to right, only a change of notation is needed to make all 
formulas apply to the case where all currents flow from right to left; and 
the case where part of j e flows to the right and part to the left can, obviously, 
occur only under conditions where the assumptions of this section are not 
fulfilled, i.e., can occur only if electrons are removed at x = or if both 
diffusion and recombination are important. For, if diffusion is negligible, 
the existence of a potential maximum at x = implies a convergence of 
electrons from both sides onto the plane x = 0, and recombination alone 
cannot annihilate electrons at a finite rate in an infinitesimal volume. 

Mention has already been made of the fact that equations such as (12) 
and (13) give an infinite density of holes when j e /j a = unfile , and are non- 
sensical for larger values of j e j a . It is easy to see why any theory which 
neglects diffusion must break down for values of j e /j a of this size and larger 
if no electrons are removed by the injection electrode. If j c /j a is too large, 
any positive field just to the right of the injection plane x = will cause 
more electrons to flow in the negative .v-direction than can be carried off 
by the current j a which flows in the region of negative x. This difficulty 
cannot be eliminated by making the field smaller in the region of small 
positive x, since making the field smaller requires a higher density of holes 
to carry the hole current j c ; and this in turn requires a higher density of 
electrons to preserve electrical neutrality. Thus, though it may be possible 
to realize experimental conditions under which the approximations of this 
section are valid for moderate values of j c /j a , increase of j e /ja above the 
critical value will always result in the building up of an enormously high 
density of holes and electrons near x = 0, and one must then consider 
diffusive transport and possibly other phenomena such as breakdown of the 
assumption that no electrons are removed by the injection electrode. 

It will be shown below that the effect of recombination on the curves of 
ni, against x at various times / can be taken into account by using a geometri- 
cal construction similar to that of Fig. 2 except that, instead of moving the 
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various points of the curve horizontally to the right with increasing time, 
one must move them along a family of decreasing curves (cf. Figs. 4, 5, 
and 6). The effect of diffusion can be described roughly as a migration of 
each point from one of these curves to another. 

3. Complete Differential Equations of the Problem 
As was mentioned in Section 1, the transport of electrons and holes along 
a narrow filament can be described by one-dimensional equations even if 
recombination at the surface of the filament causes the distribution of 
electrons and holes to be non-uniform over its cross-section. In the equations 
to follow, tih and n e will be understood to refer to averages, over the cross- 
section, of the hole and electron densities, respectively; the electrostatic 
field E can always be assumed uniform over the cross-section of the filament, 
if the latter is thin. The as yet uncertain influence of the surface on the rate 
of recombination of electrons and holes can be allowed for by writing the 
recombination rate as n R(n h /n )/T particles per unit volume per unit time, 
where R is a function which is asymptotically n h /n as its argument — K), 
and where t is the recombination time for small hole densities. For pure 
volume recombination, R = nnii e /n\ = (n h /n )(l + n h /m), while a con- 
ceivable extreme of surface recombination would be R = nh/n . 

Using this function, the continuity equations for electrons and holes can 
then be written, with inclusion of recombination and diffusion terms 

** = - • (£„„„») - S R (l) + » ( Dh p) (17) 

at dx t \m/ Sx\ ax I 

£ = ~ C4.<0 - - R (-) + g- (d. £) (18) 

dt dx t \«o/ dx \ dx ) 

where the D's are the diffusion constants, related to the mobilities /x by the 
Einstein relation 

D/n = kT/e (19) 

Using the neutrality condition n e = n + n h , subtracting (17) from (18) 
and integrating gives the equation of constancy of current, the generali- 
zation of (3) : 

£[Gu, + »)* + /*.*,] + ^ o*. - w) f; = ;(/)/*. (20) 

Solving for E gives 

j - kT(ji. - tin) -r- 
E - if Z S -u 1 (21) 
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which can be substituted into (17) to give a differential equation for tin alone: 



diih _ j 



dt 



I <^ [~ Uhllh 1 _ Ha „ />h\ 

" dx L(m« + M/.)»A + Mr»oJ T \llj 



(n + 2n,,) ~ 

, Kl 6 dx 

- — n h \x e — j r j 

e dx _{n e + n h ))i h + n e n 



(22) 



The first term on the right represents drift, the second recombination, and 
the third diffusion. This holds whether j is constant in time or not. How- 
ever, as the remainder of this memorandum will be devoted to the case 
where the currents involved are held constant after their initiation, it will 
be convenient to simplify the notation by introducing a current-dependent 
scale for * and writing the equation in terms of the dimensionless variables 



V = tlh/ttQ , S = t/r, £ = x/EqUhT = XeiloHe/jUkT 

In terms of these (22) becomes simply 



di> 
ds 



d 



D 



+ (l + M/.AO" 



i- 



«W 



+ 



(J)' 



(1 + 2.) | 



(23) 



(24) 



d£ „1 + (1+MA//0". 

where R(i>) = v(l -f v) for pure volume recombination, or = v for a surface 
recombination uninfluenced by the electron density, and where 

J = (kTe ntnl/nhT) 112 

i/o (25) 

= cr (kT/e hut) 

Numerically the characteristic field is, at 300°K, with n h = 1700 cm 2 /v sec, 3 

(kT/en, t T) m = 3.90 (t/Ims) -1 ' 2 volts/cm (26) 

Note that the importance of the diffusion term in (24) goes down in- 
versely as the square of the current density used and inversely as the square 
of the recombination time; this is because an increase in the distance the 
holes travel decreases the distance they diffuse by decreasing the concen- 
tration gradient, and also makes a given diffusion distance less serious by 
comparison with the total distance traveled. Note also that, if /i 6 = n h , 

the last term of (24) reduces simply to ( - 1 — \, but that, if /x e ^ y. h , the 

diffusion term is not a simple second derivative. 

3 G. L. Pearson, paper Q9 presented at the Washington Meeting of the American 
Physical Society, April 29, 1949. 
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4. Solution Including Recombination but Neglecting Diffusion 

It is plausible to expect by analogy with Fig. 2 that (24) can be solved, 
neglecting the last term, by a similar construction in which the curve of n h 
against x at time / is derived from that at time by moving each point to 
the right along a descending curve, instead of along a horizontal line as be- 
fore. To show that this is indeed the case, and at the same time to show 
that the diffusion term cannot so easily be taken into account, let (24) be 
written, omitting its last term, as 

* = -«(,) % - KM 
ds d% 

where <l> is just the translation into dimensionless variables of the velocity V 
encountered in (4). This can be converted into a differential equation for £ 
by writing 

(e A \ 

Wt _ "(f). 
and multiplying through by ( — J : 

©.-*>©. + •« (27) 

f dv 
orwithz* = j ^y, 

S - * = *(«0 

ds dw 

ds dw 

whence the general solution is 

| = _ J $ dw + f(s + w) (28) 

where/ is an arbitrary function. If the same transformation is tried on (24) 
with the diffusion term retained, the equation corresponding to (27) has an 
additional term on the right containing a quotient of second and first deriva- 
tives of £ with respect to v, and the simple explicit solution fails. 

To apply (28) to explicit calculation, or even to visualize it physically, it 
is necessary to determine the proper form of the arbitrary function / to fit 
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the boundary conditions of the probfem. This is most conveniently done by 
introducing a family of curves as suggested by the analogy with Fig. 2. The 
analogy suggests that we should try to find curves in the v, £ plane (the full 
curves of Fig. 4) such that a point can move along any one of them with 
velocity components 

§ = *> £--*■ 

as as 

The equation of any such curve is 



dp 
or 



$ - -*/* 



«fr, *) = [' *~£ (29) 

where p , the intercept of the curve on the c-axis, is taken as a parameter 
distinguishing the curve in question from others of the family. A point which 
starts at height v on the v-axis at time s = will reach height v at time 

*,*)- (* %• (30) 

Thus, after time s, the locus of all points which start at all the various 
heights v will be the curve obtained by eliminating i> between (29) and 
(30) (shown dotted in Fig. 4). That this curve is, in fact, of the form (28) 
and therefore a solution of the differential equation is easily seen by writing 
(29) and (30) in terms of integrals taken from some arbitrary but fixed 
lower limit: 

As fn is varied both the integrals with upper limit v u will vary, and either 
can be expressed as a function of the other: 



whence 



n?-f(r d i) 



which is identical with (28). 
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The equations (29) and (30) of course apply only to the portion of the 
curve of v against £ which is derived from starting points v on the v axis 
which are less than the maximum value Vi corresponding to the value -nui 
given by (12): The points for v < v y are merely initiated at time s = T).a.nd 
propagated by the differential equation from then on; the point v =■ n , 
£ = 0, on the other hand, remains a source at all times from the initiation 
of the injection onward. Thus the complete curve of v against £ for any 
positive ^ follows the dotted construction of Fig. 4 from the £ axis up to 
where it intersects the full curve corresponding to v = n , after which it 




^ccx 



Fig. 4— Schematic illustration of the method of constructing the curve of hole density 
n h against distance x from the injection electrode at some given time, taking account of 
recombination but neglecting diffusion. 

follows the latter curve, as indicated by the crosses in the figure. The steady- 
state distribution is thus simply the full curve for v = v x . 

For explicit calculation for the case of pure volume recombination one 
must insert $ = 1/[1 + (1 + Ph/M\ R = "(1 + ") into (29) and (30). 
The integrations are easily carried out and give 

i = [ t , 1 +." ,/ ^ ; + In 



n 



+ (1 + Ma/m.) 



^o 



1 4- (1 + Ph/p.) 



v« 



+ *-, 

Ma 2 



In 



1 + (1 + 



1 + vo 



Hh/fte) vo I 
VO J 



— [same with v instead of vo] 



= In 



VQ 



1 + Vo 



- In 



1 + v 



(31a) 
(32a) 
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For the case of a surface recombination uninfluenced by electron concen- 
tration one obtains similarly, with R = v: . ■ 

* = |_1 + (1 + hMh " ^ 1 + (1 + WlOfJ 

— [same with v instead of v Q ] (31b) 

s = In -° (32b) 



When /x e = 3/X/./2, as for germanium, (31a) and (31b) become respectively 



, V 2 , , *> , 9 . 1 + 5i*/3 

£ = 1 i c /i + ln i ■ c it + 7 In 



h 5i*/3 1 + 5p /3 '4 1 + vq J 

same with v instead of vq] (33 a) 



and 



€ „,-_ +ln ' 



- r— i 

|_1 + 5v /3 ' 1 + 5po/3. 
— [same with j» instead of i/ ) (33b) 



These can also be written, using (32a) and (32b), 

, _ , 5/2 5/2 5 [ (I + 5,0/3X1 + y ) l 

* f l + 5vo/3 " 1 + Sr/3 + 4 ^ L(l + 5p/3)(1 + v )j C34a ' 

and 

/f + 3/S \ 1 1 

* = 5 + ln U+vsJ " rrv3 + r+w3 (34b) 

Figures 5a and 5b show as a full curve the plot of eq. (33a) for the case 
v = oo, and the full curve in Fig. 6 shows in the same way the plot of 
(33b) for i> = co. Changing v of course merely shifts either curve hori- 
zontally. Note the very sharp increase of v for small £, which shows up in 
pronounced manner on the expanded scale of Fig. 5b. The corresponding 
values of s, computed from (32a) or (34a), are marked on the curve of Fig. 5; 
the corresponding marks on the curve of Fig. 6 also represent values of 5 at 
intervals of 0.2, but are not labeled with absolute values because (32b) is 
infinite for i»o = co . 

For large £, v becomes very small and it becomes legitimate to expand the 
logarithms. The first few terms of the resulting asymptotic expression for £ 
are, for v = co and the recombination function leading to (31a), 

£ ~ (g - lj In (1 + jia/m.) - (1 + nM ~ In v 

+ (3 + 2 M a//i.> (35a) 
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Fig. 5— Steady-state curve of hole density m, against distance .v, for the case of ideal 
volume recombination (recombination rate = »i»./T»o), and asymptotic approximations 
to this curve. 



EXCESS SEMICONDUCTOR HOLE TRANSPORT 417 

while, for the recombination function leading to (31b), 

£ „ _ i n(1 4- M ,,/ Me ) -l-]w+ 2(1 4- nMv (35b) 

In Figs. 5a and 6 the lower dotted curve represents the sum of the terms of 
(35a) or (35b) respectively as far as the term in In v: in this approximation 
the dependence of v on £ is exponential. An exponential behavior of this 
sort is assumed in the small-signal theory of the modulation of the resistance 
of a filament of semiconductor by hole injection.' 1 The upper dotted curve 
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Fig. 6 — Steady-stale curve of hole density tih against distance x, for the case of ideal 
surface recombination (recombination rate = «a/t), and asymptotic approximations to 
this curve. 

in Figs. 5a and 6 is a plot of (35a) or (35b), respectively, with the linear 
term included. It will be seen that in both figures the simple exponential 
approximation is already quite far off when v = iih/Ho = 0.1, though it 
improves rapidly for smaller v. 

Figure 7 shows a sample plot of v against £ for the case of ideal volume 
recombination (eqs. (31a) etc.), for the numerical conditions s = 1, v\ «* 0.3 
(cf Fig. 4). According to (12), whose validity at £ = is unimpaired by the 
occurrence of recombination, this value of v x implies /„'/'< = 6.5. The left- 

* \Y. Shockley, G. L. Pearson, and J. R. Haynes. Bell Sys. Tech. Jour., this issue. 
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hand portion of this curve is simply traced from Fig. 5, with a horizontal 
shift sufficient to give an intercept at v = 0.3; the right-hand portion was 
constructed by placing the paper for Fig. 7 over that for Fig. 5, shifting 
horizontally until the point corresponding to one of the values of 5 marked 
on Fig. 5 lay on the axis of ordinates of Fig. 7, marking the position of the 
point labeled with one plus this value of s, and repeating. 
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Fig. 7 — Variation of hole density «/, with distance x at time I = r assuming mm = 0.3 no 
recombination rale = HUli/rno, and neglecting diffusion. 

Figure 8 shows sample plots of v against s for the same case of ideal 
volume recombination, with vi = 0.3, for £ = 0.5 and £ = 1.0. Curves for 
a different v\ would start out exactly the same, but rise higher. The rising 
portion of the curve for £ = 0.5, for example, was constructed from the 
curve of Fig. 5a by locating various points (£, v) on the latter curve and 
associating with the v value of each such point a value of 5 equal to the 
difference of the 5 values marked on the curve of Fig. 5a for the two points 
abscissae £ and (£ - 0.5). As Fig. 5a was prepared entirely by slide rule, 
the accuracy is not all that can be desired; the individual computed points 
are shown to give an idea of the magnitude of the computational errors. 
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For convenience in future calculations the equations will be appended 
which correspond to (31) to (34) when, instead of Hi, , the field E is used as 
dependent variable in the differential equations. In terms of the dimension- 
less variable 



€ = E/E = 



1 + V (l + Hh/He) 



(36) 



and the parameter e corresponding to v = v , the equations are, for ideal 
volume recombination (eqs. (31a) etc.), 



i = [(1 + m ./ma)«o - T In (l + ^ ft) + In (1 - eo)l 



— [same with e instead of eo] 



. (1-ft) . (1-0 

5 = In "7 \ — In 



( 1+ £«) 0+2-) 



(37a) 



(38a) 



while, for the recombination function leading to eqs. (31b) etc., 

1 - eo 



£ = eo — e + In 



5 = In 



1 - e 



1 


"1 




1 


eo 




1 






1 


_e 





(37b) 



(38b) 



The electrostatic potential U is 

U = — / E dx — — eIuht I e </£. 
In the steady state the relation between e and £ is given by (37) with e 

(Mc + V-h) je - ,, . 

-r— — rr . For this case 



set equal to ei which, by (13), is 1 — 



Ma Ua + je) 



U = —Eqixht «$ 



= —EoflhT 



-M 

c(m!/ma - 



1) - e 2 (l + m./ma)/2 



(39a) 



- In (1 - «)■ - t In f 1 + ~ e)| + const. 
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for ideal volume recombination; while, for the assumptions leading to eqs. 
(31b), etc. the relation is 



V = -Eluh t[« - e 2 /2 - In (1 - e)] + const. 



(39b) 



Thus, in the steady state, the difference in potential between any two points 
to the right of x = can be obtained by finding the values of e for these 
two points by (37), and then using these to evaluate the difference in the 
values of (39) at the two points. To the left of x = 0, of course, E is constant 
and equal ioj a /<r. 

5. Diffusion Effects 

Diffusion will obviously be very important at small values of £ = x/E n h r 
when n,,i is large, because of the tremendous concentration gradients which 
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Fig. 8— Transient behavior of »«, with time at position x/Eo hut = 0.5 and 1.0, as- 
suming wai = 0.3 Ho and recombination rate = »aM«/t«o. 

Figs. 5 and 6 predict for such cases. Also, of course, diffusion will round off 
the discontinuities in slope which appear at the front and rear of the transient 
disturbance as in Fig. 7 and Fig. 8. At other points the importance of diffu- 
sion effects can be roughly estimated either by comparing the diffusion cur- 
rent with the drift current or by comparing the divergences of these two 
contributions to the current, i.e., the last and first terms on the right of (24). 
Referring to these terms in (24) we have 

[diffusion currentl _ /f\ 
drift current J \j/ 

[ div. diffusion current "! 
div. drift current J 



[1 +2 



91 S] 



(40) 



-(? 



, (41) 



[l + (i + m/riwULi + 2u) €l 



TJT + (1 - Hh/fie) 77 
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For the steady-state curve approximate values of the expressions (40) and 
(41) can be computed by evaluating the derivatives of £ with respect to v 
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Fig. 9 — Asymptotic magnitude of diffusion terms in the steady-state flux of holes, 
when jl J is large. 

from (29) or (31). For the case of ideal volume recombination with m^/ma = 
3/2 this gives, if the diffusion effects are not too large, 



piiffusion currentl ^ /VY , 

|_ drift current J \ / / 

(pj (1 + *f)'(1 4- V* + 2U 2 4- -W). 



div. diffusion current"! 
L div. drift current J 



(43) 



These functions are plotted in Fig. 9. From this figure one can estimate 
roughly when diffusion will begin to have serious effects other than a slight 
rounding of the leading and trailing ends of the transient. For example, if 
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it is desired that the ratio (43) be less than about 0.1 in the steady state 
for values of v as high as 0.3, the upper curve of Fig. 9 shows that the current 

density used must be large enough to make f - J < t~, i.e., j > U.7 J, 

where J is given by (25) and (26). 

An approximate evaluation of (40) and (41) in the transient region can 
be performed by graphical or numerical differentiation of a curve such as 
that of Fig. 7. For example, a rough calculation based on Fig. 7 gives, in 
the middle of the transient portion (£ = 0.75), 



[ div. diffusion current "! ^ ~ /i V 
div. drift current J" \jj 



More important and also more difficult to estimate is the effect of diffusion 
in rounding off the front and rear edges of the transient. Various ways can 
be devised to estimate a rough upper limit to the amount of rounding off 
to be expected. One such is to compute what the diffusive flux just behind 
the front of the advancing disturbance would be if the distribution of holes 
were the same as in the absence of diffusion. Under conditions where diffusion 
is not too serious the time integral of this diffusive flux between any two 
times can be equated to the increase in rounding of the front, as measured 
by the area between an ideal curve such as that of Fig. 7 and the actual 
curve of v against £ for the same time s. The integration cannot be extended 
back to time zero, however, since the integral of the flux diverges logarith- 
mically. The fact that the diffusive flux is actually finite instead of infinite 
of course arises from the fact that at small times the concentration gradient 
a short distance behind the front can no longer be approximated by the 
gradient which would obtain in the absence of diffusion, but instead is very 
much less. Thissuggests that an upper limit to the total diffusive flux passing 
into the region of the front from time to time 5 can be obtained by taking 
the flux computed as described above between the times s and s, and adding 
to it the total number of holes which have left the injection electrode be- 
tween time and time s . Since this gives an upper limit for any s , one 
may use the minimum of the resulting sum as s is varied. 

The results of some sample calculations of this sort are shown in Fig. 10, 
which refers to the same time, currents, and recombination function as Fig. 7, 
viz., s = I/t = 1.0, j e /j a = 2/13, ideal volume recombination. The full 
curve is the transient portion of Fig. 7 replotted on a larger scale. The lower 
dotted curve is a curve drawn in by hand in such a way as to make the 
area between it and the full curve equal the upper limit computed in the 
manner just described, for the case j = 1007. The upper curve was drawn 
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similarly for 7 = 31.6/. Since the true curve of v against £ must lie between 
the dotted curve and the full curve in each case, it can be concluded that 
for times and current ratios of this order the diflusionless theory of Section 4 
gives a useful approximation to the transient when j > 1007. At the other 
end, it seems likely that iorj < 10 J the theory of Section 4 has no quantita- 
tive utility at all in the transient region. 

When diffusion effects are sufficiently great, account must also be taken 
of the fact that the boundary conditions at the injection electrode (x = 0) 
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Fig. 10 — Approximate magnitudes of the rounding of the front by diffusion for various 
values of j/J, for the case/ = r,j,/j a = 2/13, ideal volume recombination. Ordinate is 
proportional to hole density, abscissa to distance from injection electrode. 

take a different form from those in the absence of diffusion. In the absence 
of diffusion and with the assumption that only holes are injected at x = 0, 
the current just to the right of x = must consist of a contribution j e from 
holes and a contribution j a from electrons, while the current just to the left 
of x = is purely electronic and of magnitude j a . This implies, as we have 
seen in Section 2, that the hole density be discontinuous, with the value ni.\ 
given by (12) just to the right of * = 0, and the value zero just to the left. 
But if diffusion is allowed, the hole density must be continuous. For the 
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idealized case where holes are injected on the plane x = and no electrons 
are removed there, the equations to be satisfied are 

«.(S) + + «.^- A (£). + -**- (44) 

-Ih (^') + + n htth E + = j - c - D h (%£) + *«*. (45) 

De (j£) ~ Dh (j^) +<** + *»)*-*'£/« (46) 

where subscripts + and — refer to conditions just to the right of X = 
and just to the left, respectively. Using the neutrality condition n e = n + «* 

these are three equations for the five unknowns f -~ J ,E ± , n h . To complete 

the determination of these quantities the differential equation (22) must be 
solved and the boundary conditions imposed that n h — * as x — * ± <». 
Actually the problem of estimating conditions at .v = may not be quite 
as formidable as the preceding paragraph suggests, at least if the diffusion 
parameter J/j is reasonably small and \ij e /j a is also not too large. For such 
cases the "upstream diffusion" of holes into the region of negative x will 
probably reach a steady state in a very short time. Solutions of the steady 
state differential equation in this region have been obtained numerically by 
W. van Roosbroeck (unpublished). Such solutions will give one relation be- 
tween ;//, and f — - ) ; another relation, in the form of a fairly narrow range 

of limits, is provided by the fact that f — ? 1 will under these conditions be 

« ( "l J f being in fact probably somewhere between zero and the value 

for the diffusionless case with the same value of uia . 

Of course, if the mathematical solution for this one-dimensional idealiza- 
tion is to be applied to a case where holes are injected into a filament by a 
pointed electrode on its boundary, little meaning can be attached to vari- 
ations in the «/, of the mathematical solution within a range of x values 
smaller than the diameter of the filament. 

6. Summary and Discussion' 

There are three principal factors which limit the range of conditions 
within which the present theory provides a useful approximation to the 
transient behavior of w* as a function of / and x. These are diffusion, trap- 
ping, and departure from one-dimensional geometry. If the geometry is 
sufficiently nearly one-dimensional and trapping is negligible, the discussion 
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of Section 5 shows that the theory of Section 4, with its neglect of diffusion, 
will give a useful approximation to the truth whenever the field in which 
the holes migrate is sufficiently strong — e.g., strong enough to make the 
current density /* > 100 ./, where J is given by (25) and (26). The obtaining 
of "sufficiently strong" fields without excessive heating or other undesirable 
effects is facilitated by the use of specimens with as long a recombination 
time r as possible, and by the use of specimens of low conductivity. How- 
ever, it is hard to say how low the conductivity can be made without danger 
that the "no trapping" assumption will break down, since for this assumption 
to be valid the density of hole traps must be « the density of donors. 

The numerical predictions of the theory depend upon the way in which 
the rate of recombination is assumed to depend upon the concentrations of 
electrons and holes, i.e., upon the form of the function R(v) introduced in 
(17) and (18). The full curves of Figs. 5 and 6 give the steady-state depend- 
ence of Uh on x for two simple assumptions regarding R(u), the dependence 
corresponding to any given boundary value nu at x = being simply ob- 
tained by a suitable horizontal shift of the curve plotted. When the currents 
are held constant after their initiation, the auxiliary time scale in these 
figures can be used to construct the transient disturbance at any time, by 
the methods described in connection with the examples of Figs. 7 and 8. 

These results should hold for a plane-parallel arrangement of electrodes or, 
to a good approximation, for electrodes placed along the length of a narrow 
filament, provided the »;, appearing in the equations is interpreted as a cross- 
sectional average of the hole density and provided the other assumptions 
given in Section 1 are fulfilled. It is easy to see, however, that practically 
the same equations apply to cases of cylindrical or spherical geometry, in 
the approximation where diffusion is neglected. For, in these cases, the 

d 1 a 

original equations (17) and (18) merely have ~ ( • • • ) replaced by ~~ (r ■ • ■ ) 

ox r or 

1 o ... . 

or ~l T~ (r- • • •); if the diffusion terms are neglected the solution is the same 
T or 

as before with x replaced by r 2 /2 (cylindrical case) or r*/3 (spherical case) 
and with j replaced by I/2ird (cylindrical case, d = thickness of sample, 
/ = total current) or by I /Air (spherical case). However, it may be difficult 
to realize experimentally conditions approximating cylindrical or spherical 
geometry which satisfy the requirement that diffusion effects be small. 

Another generalization which is easily made is the removal of the assump- 
tion that no electrons are withdrawn by the electrode at x = 0. As far as 
conditions to the right of .v = are concerned (Fig. 1), the only change 
required in the diffusionless theory is to interpret j e as the current density 
leaving the emitter electrode in the form of holes, rather than as the total 
current from the emitter electrode, and to interpret j a as the sum of the 
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current leaving the emitter electrode in the form of electrons and any current 
to the left of * «= 0. 

It should also be clear that the entire analysis of this paper, though it has 
for definiteness been formulated for the case where holes are injected into an 
excess semiconductor, applies just as well to any case where electrons can 
be injected into a defect semiconductor. For the latter case it is merely 
necessary to interchange the subscripts e and h in the formulas. Though the 
types of experiments discussed in this paper have to date only been reported 
for w-type germanium, the occurrence of similar phenomena in />-type speci- 
mens is indicated by the successful use of such specimens in transistors. 6 

An interesting and possibly quite useful phenomenon should occur when, 
after establishment of a steady state, the current j e is suddenly turned off. 
There will result a transient disturbance propagated in the direction of in- 



(a) 




Fig. 11— Schematic variation of hole density m, with distance x, illustrating formation 
of a shock wave by quickly decreasing j c to zero, for the case where j = j, + j a is kept 
constant. 

(a) Immediately after reduction of j, to zero. 

(b) Later time. 

creasing x, which is very much like a shock wave in a gas. This, the most 
interesting feature of the phenomenon, will occur regardless of whether j a 
remains constant when j, is cut off; however, the simplest example for il- 
lustrative purposes is the case where j a is increased by the amount^ at the 
instant when the latter is cut off, so thaty remains constant. For this case, 
illustrated in Fig. 11, the values of n h ahead of the advancing front will 
remain the same at each point as in the previous steady state. Just behind 
the front, n h must drop abruptly to zero. If j/J is large, where / is given 
by (25), the drop will be extremely sharp. For the change in the form of the 
front with time is compounded out of diffusion and propagation with variable 
velocity along descending curves, as shown schematically in Fig. 4. Since 
the latter propagation involves a more rapid motion to the right, the smaller 
iih , it tends to steepen the front, and this steepening must continue until 

6 W. G. Pfann and J. H. ScafT, paper presented at the Cambridge Meeting of the Ameri- 
can Physical Society, June 16-18, 1949. 
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the diffusive spreading becomes sufficient to counterbalance it. It is not 
necessary, for the production of a steep front of this kind, that the decrease 
of j e to zero be brought about with corresponding rapidity; even a gradual 
decrease of j e will lead to a front which becomes steeper as it advances, 
and if the decrease of j e is not too gradual a "shock front" will have devel- 
oped after a short distance. The order of magnitude of the "shock front 
thickness" can be estimated by finding the value of the time At for which 
the diffusion distance Axd = (2D A /) equals the difference Axy between 
the drift distances of the holes at the top and bottom of the front, i.e., 
Ax r = [V(0) — V()h)]Al, where V is given by (4) and tth is the height of 
the front. For this value of At, 

Ax D = 2D/[V(0) - V(n lt )] (48) 

and this is presumably of the order of magnitude of the thickness of the 
front. If D is interpreted as Dn = kTn h /e, which is good enough for the 
present purpose, this gives 

2kT 1 

A.v D = ^t- — -T,. (49) 



Ti- i T 

L 1 + v(l + Hh/vJI 



Of course, this extremely sharp front can be realized only when the condi- 
tions of one-dimensional geometry are accurately fulfilled. When the geome- 
try is made sufficiently ideal, observation of the thickness of the "shock 
front" can provide a valuable check on the validity of the basic assumptions 
of the theory such as the neglect of trapping. 6 

The author would like to express his indebtedness to many of his col- 
leagues, and especially to J. Bardeen, J. R. Haynes, and VV. van Roosbroeck, 
for many illuminating discussions of the topics covered in this paper. 

6 The accompanying paper by W. Shockley, G. L. Pearson and J. R. Haynes describes 
some observations of this shock wive effect, though under conditions where v «C 1, so 
that the thickness of the front as given by (49) is still fairly large. 
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